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Abstract 

In this paper we study the experimental consequences of a the- 
ory which naturally has a heavy, stable (or almost stable) gluino. 
We define the boundary conditions at a messenger scale M ~ 10 14 
GeV which lead to this alternative phenomenologically acceptable 
version of the minimal supersymmetric standard model. 

In this theory, either the gluino or the gravitino is the lightest 
supersymmetric particle [LSP]. If the gravitino is the LSP, then 
the gluino is the next-to-LSP with a lifetime on the order of 100 
years. In either case, the gluino is (for all practical purposes) a 
stable particle with respect to experiments at high energy accel- 
erators. Thus the standard missing energy signature for SUSY 
fails. A stable gluino forms a color singlet hadron, the light- 
est of which is assumed to be an isoscalar gluino-gluon bound 



state (Rq). The Rq has strong interactions and will interact in a 
hadronic calorimeter; depositing some fraction of its kinetic en- 
ergy. Finally, in the case the gravitino is the LSP, bounds from 
searches for stable heavy isotopes of hydrogen or oxygen do not 
apply to the metastable Rq. 
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1 Introduction 



Supersymmetry [SUSY] is a strongly motivated candidate for new physics 
beyond the Standard Model [SM] . The minimal supersymmetric particle con- 
tent is well defined and the interactions of all the new superparticles [sparti- 
cles] are constrained by the observed SM interactions as long as the theory 
has a conserved R parity. In order to search for SUSY one necessarily focuses 
on specific signatures which may rise above the standard model backgrounds. 
These signatures depend on how supersymmetry is broken; i.e. on the spar- 
ticle spectrum. There are two classes of models for SUSY breaking which 
have been studied extensively. They are the minimal supergravity model [TJ 
(or constrained minimal supersymmetric model [CMSSM]0) defined by 5 
parameters and low energy gauge-mediated SUSY breaking [GMSB]|| |J 
defined by 3 parameters. The 5 parameters in the CMSSM model are a 
universal scalar mass mo, a universal gaugino mass M1/2, a trilinear scalar 
coupling proportional to A and /i, fiB = m§ in the Higgs sector. In this 
theory, the LSP is the lightest neutralino or a sneutrino. The dominant sig- 
nature for SUSY is missing energy due to the escape of the LSP from the 
detector. In low energy GMSB, the effective scale of SUSY breaking is given 
by A = F/M 10 5 GeV. The other two parameters are /x, fiB = m\ 

as before. In this case, the gravitino with mass of order 10 - 1000 eV is the 
LSP. The next to lightest SUSY particle [NLSP], for example the lightest 
neutralino, can decay into a gravitino and a photon. This can lead to en- 
ergetic photons plus missing energy, such as the single e + e~77 + missing 
energy event seen at CDF. 

In a recent paper || it was shown how GMSB with Higgs-messenger mix- 
ing in an SO(10) theory naturally leads to a gluino LSP. Since the gluino 
carries color, it will be confined. The stable LSP is then assumed to be a 
gluino-gluon bound state, a glueballino {Rq). Rq is a hadron and thus in- 
teracts strongly in a hadronic calorimeter. As a result the standard missing 
energy signature for SUSY fails. This theory has a characteristically different 
spectrum of squark, slepton and gaugino masses. In this paper we analyze 
the low energy spectrum of the theory, starting with boundary conditions at 
the messenger scale M ~ 10 14 GeV, using two (one) loop renormalization 
group equations for dimensionless (dimensionful) parameters down to the 
weak scale, mz- We self-consistently require electroweak symmetry break- 
ing using the one-loop improved Higgs potential. The model has six SUSY 
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breaking parameters — A ~ 10 5 GeV, a ~ b ~ 1CT 1 — 1CT 2 , //, /iB = m| 
and d (the size of the D term contribution in units of M2). One of these 
parameters, b, allows us to arbitrarily vary the gluino mass. 

Finally we note that the version of the model presented in this paper 
includes the SUSY breaking contribution of an anomalous U(l)x, in addition 
to GMSB at the messenger scale M0. We show that the additional D term 
contribution is necessary in order to obtain a phenomenologically acceptable 
theory. In particular, without the D term the light Higgs boson h would 
be unacceptably light. In the following, we assume that the contribution 
to scalar masses from GMSB and the U(l)x D term are comparable. In a 
recent paper ||, we have presented a model which dynamically breaks SUSY 
and leads to comparable SUSY breaking effects from gauge-mediated and D 
term sources. The messenger scale in this theory is determined by the Fayet- 
Iliopoulos D termpj of U(l)x- 1 This theory suggests that the boundary 
conditions we consider may be "natural." 



2 Boundary conditions at the Messenger scale 

The boundary conditions at the messenger scale are determined by two 
sources of SUSY breaking — gauge-mediated and D term. 



2.1 Gauge mediated SUSY breaking 

The messenger sector is defined by an SO(10) invariant superspace potential 
W D A a 10 H A 10 + X s S 10 2 

1 The magnitude of the Fayet-Iliopoulos D term £ is given by £ = e M Pl where e = 

9 *92 < k 2 ( g^Mp; ) 2 an< ^ * ne ^ring scale M st is determined by the compactification scale. 
In weak coupling M st — 22L Mp\. In the strong coupling limit of the E s x Eg heterotic 
string^], on the other hand, the compactification scale is identified with the GUT scale. 
Note, however, in a recent paper j^] it has been argued that even in the strong coupling 
limit of the heterotic string, the scale of the D term is set by the Planck scale and not the 
compactification scale. This argument has some caveats as discussed in the paper. For 
example, it might not apply if the anomalous gauge symmetry eminates from D brane 
modes. Another possibility is that the anomalous U(l) in the effective four dimensional 
theory is a linear combination of U(l)s on the two 10 dimensional boundaries of the E$ x Eg, 
heterotic string. In this case the axion which cancels the anomally is partially in the S and 
T moduli of the theory. As a consequence the scale of the anomally becomes arbitrary. 
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+Ai 16i A 16i + A 2 16 2 A 16 2 + A S 16i 16 2 (1) 

The adjoint A is assumed to get a vacuum expectation value [vev] 

< A > =(B-L) M G (2) 

where B — L (baryon number minus lepton number) is non- vanishing on color 
triplets and zero on weak doublets in the 10 and the singlet S is assumed to 
get a vev 

< S(x, 0)> =S + 6 2 F s . (3) 

We take S ~ 10 14 GeV with the effective SUSY breaking scale in the observ- 
able sector given by 

A = F s /S ~ 10 5 GeV. (4) 

The field 10h includes the weak doublet and color triplet Higgs, while 10 
contains the minimal set of messengers. The first two terms in W correspond 
to the natural doublet-triplet splitting mechanism in SO(10)HlQ|. Note that 
the doublet messengers obtain mass at the scale S, while the triplet mes- 
sengers obtain mass at the GUT scale. This theory also includes a natural 



mechanism for suppressing proton decay when S « Mg WM - Finally, the 
last three terms in the superspace potential include four additional sets of 
messengers. They are introduced solely to give gluinos mass at one loop.0 
Note, the R symmetry which keeps gauginos massless is spontaneously bro- 
ken at S. 

In this theory, gauginos obtain mass at one loop 

= sum b 2 a 

y tt 
M 2 =^g^A(l + 46 2 ) 

Mi =|^A(l + f6 2 ) (5) 



where the parameter b = A S/ vAi A 2 Mq is derived from the last 3 terms in 
W. The parameter b is naturally of order 1CT 2 A/VAi A 2 . For the purposes 

2 The first two terms do not contribute to gaugino masses at one loop due to an acci- 
dental cancellation as shown by one of us (K.T.). 
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of our analysis, we take b = 0.1. However b is clearly a free parameter which 
may be used to vary the gluino mass. 
Scalars obtain mass at two loops 

m 2 =2A 2 {C 3 (^)V + 4& 2 ) 

+C 2 (^) 2 (l+46 2 ) 

+^i(^) 2 (I + I« 2 + 46 2 )} (6) 

where C 3 = 4/3 for color triplets and zero otherwise, C 2 = 3/4 for weak 
doublets and zero otherwise, G x = §(V/2) 2 and a = X s S/X a M G . Note, the 
parameter a derives from the first two terms in W. In our analysis, we take 
a = 0.01. The arbitrary parameters a and b are clearly independent. 
Finally the gravitino mass is generically given by 

m 3/2 =F s /V3M Pl (7) 

with the reduced Planck scale, M Pl = l/y/8nG N = 2.4 x 10 18 GeV. Note, 
supergravity corrections to squark and slepton masses are proportional to 
the gravitino mass. 

In the particular model of SUSY breaking discussed in ||, the field which 
gets both a scalar and F component vev is the third component of an SU(2) p 
vector field, S 3 . In this theory, the field S is a composite superfield, with 
S = S 3 2 /M st and F s = 2S 3 F S3 /M st . The scale M st is associated with a 
string scale of order M G - Hence with a maximum value of S3 ~ 10 15 GeV we 
obtain S ~ 10~ 2 M G ~ 10 14 GeV. The upper bound on S 3 ~ 10 15 GeV has 
been assumed in order to suppress flavor changing neutral current processes 
induced by supergravity corrections to squark and slepton masses. With the 
gravitino mass now given byf\ 

m 3/2 =F s jV3M Ph (8) 

the ratio of the supergravity contribution to squark and slepton masses (fixed 
by the gravitino mass) to the GMSB contribution scales as 

m 3/2 /M 2 =-^-(S 3 /M Pl )~W 2 (S 3 /M Pl ) -0.04. (9) 

3 The gravitino mass is set by the largest SUSY breaking vev in the theory. The D 
term contribution is of order the weak scale and thus is negligible. In the particular SUSY 
breaking sector of the theory given in || , the fundamental SUSY breaking vev is given by 
Fs 3 — Yst ~ 10 Fs- This explains the change in going from cqn. ^ to ||. It is clear 
that the value of the gravitino mass is model dependent. 



6 



Hence with S3 < 10 15 GeV, supergravity gives at most a 0.16% correction to 
squark and slepton masses squared. This may be sufficient to suppress large 
flavor changing neutral current processes. 

As an example of the spectrum at M we have (for A = 10 5 GeV) 

m g = 13 [b/.lf GeV 
M 2 = 330 GeV 

m 3/2 = 12 (S 3 /10 15 GeV) GeV. (10) 

In addition, scalar masses are fixed by the largest SUSY breaking scale in 
the problem, M 2 ; with right-handed squarks and sleptons being the lightest 
scalars obtaining mass only via weak hypercharge interactions. 

Thus, in this model the gluino (or possibly the gravitino) is naturally the 
LSP. The gluino mass is suppressed by two powers of S/Mq- One power 
comes because the color triplet messengers have mass of order Mq and the 
second comes because the R symmetry in this sector of the theory is only 
broken at the lower messenger scale S. 

2.2 D term SUSY breaking 

For phenomenological reasons we assume that SUSY is also broken by the D 
term of an anomalous U(l)x gauge symmetry. In a recent paper|| we have 
shown that it is possible to obtain D term SUSY breaking and GMSB with 
contributions to scalar masses of the same order. In this paper we simply 
parametrize the D term contribution to scalar masses by the formula 

5 D m 2 a =Q x a dM 2 2 (11) 

where is the U(l)x charge of the field labelled by the index a and d is 
an arbitrary parameter of order 1 which measures the relative strength of D 
term vs. gauge mediated SUSY breaking. 

Since we are working in the context of an SO (10) GUT, Q x necessarily 
commutes with the SO(10) generators. In order not to consider a random 
U(l) gauge symmetry we shall assume that U(l)x is embedded into E 6 with 

Qf =1 -2 4 

on 16 10 1 c27ofE 6 . (12) 
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This U(l)x is naturally family independent and is thus safe from inducing 
large flavor changing neutral current processes. It is also well-motivated in 
the context of string theories [[T^]. 



2.3 Summary of boundary conditions at M 

Thus the boundary conditions at the messenger scale M from GMSB (for 
scalar and gaugino masses) and from D term SUSY breaking (for scalar 
masses only) are given in eqns. |5|, |], and [11. There are 3 additional param- 
eters — the supersymmetric Higgs mass /i, and the two soft SUSY breaking 
parameters, i.e. a scalar Higgs mass m^ 2 and the trilinear scalar coupling A. 
To leading order we have 113 1 



A =0 . (13) 

On the other hand, for the purpose of this analysis, we leave \x and m^ 2 arbi- 
trary. In principle, these parameters must be determined by new physics 
which we have not considered in this analysis. Using the 6 parameters 
A, a, b, d, /i, m 2 defining the boundary conditions at M — 10 14 GeV 
(eqns. ||, |6], |TT] and [13]; /i and m 3 2 arbitrary) we renormalize the full set of 
dimensionful (dimensionless) parameters at one (two) loop down to Mz- /i 
and m\ are fixed by requiring electroweak symmetry breaking. For the bot- 
tom mass, we have included three loop QCD running of m&(/x) from /i = mj 
to mz and one loop SUSY threshold corrections at mz, relevant for large 
tan MR . 



3 Sparticle masses at the Z scale 

3.1 Why introduce D term SUSY breaking? 

Consider first the low energy spectrum in the absence of D term scalar masses, 
i.e. d — 0. In fig. 1, we show the squark, slepton, Higgs and gaugino masses 
for the case M 2 = 110 GeV, \i < 0, a = 0.01, and b = 0.1 plotted as a 
function of tan/3. Note, the one-loop corrected light Higgs mass (the solid 
brown line) is below the present experimental lower bound, m^ > 90 GeVjUj. 
We find mh < 60 GeV. This is due to the fact that at tree level the Higgs 
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200.0 



M 2 = 110 GeV, b = 0.1, (i< 



150.0 



100.0 



50.0 




0.0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 1 

1.50 1.60 1.70 1.80 1.90 2.00 

tanp 

Figure 1: Mass spectrum of squarks, sleptons, neutralinos, charginos, gluino, 
and Higgs for the case M 2 = 110 GeV, \i < 0, a = 0.01, and b = 0.1 plotted 
as a function of tan/?. The right-handed up squarks (green), down squarks 
(red), sleptons (pink), neutralinos (blue), charginos (black), and gluino (dot- 
dashed green) masses are given by the respective colored lines. The first and 
second (third) generation squark and slepton masses are represented by solid 
(dashed) lines. The lightest Higgs mass is given by the solid (dashed) brown 
line at one loop level (tree level). 
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mass satisfies the approximate mass formula 

m h ~ m z cos 2/3 (14) 

with tan/3 < 2 in fig. 1. We could in principle increase the Higgs mass by 
increasing tan /3, however, due to the relation 

^ - 2~ H tan^/3-1 

as tan /3 increases, /z decreases and as a consequence the chargino mass also 
decreases. The lower bound on the chargino mass > 85 GeV[|l(| is 
violated for tan/3 > 1.83.0 Of course, we could in principle increase the 
chargino mass by increasing M 2 , however, the right-handed SUSY breaking 
stop mass squared is negative at the weak scale and is driven further negative 
when M 2 increases since this is the dominant source of SUSY breaking. For 
M 2 > 150 GeV, the lightest stop mass becomes less than (see fig. 2). 

Thus we cannot make M 2 very large. Hence since both tan/3 and M 2 are 
severely constrained from above, the light Higgs mass is always below the 
experimental bound. We have checked that the problem of a light Higgs, 
chargino and stop is even worse for /i > 0. 

In order to solve this problem we have been forced to consider additional 
soft SUSY breaking mass contributions. In order not to lose predictability, we 
have considered the well motivated addition of an anomalous U(l)x and its 
associated D term. We now consider how this one new contribution resolves 
the light Higgs problem. 

3.2 The addition of D term scalar mass corrections 

As discussed, we parametrize the D term contribution to soft SUSY breaking 
scalar masses at the messenger scale by 

5 D ml =Q*dM 2 2 (16) 

with 

Q x = 1 (17) 

4 This bound assumes the standard missing energy signature for SUSY. This limit is 
not applicable in our case. However, there will be a lower bound on the chargino mass 
coming from the visible Z width or e + e _ — > hadrons which is comparable. We will discuss 
these limits on the chargino mass within our framework shortly. 
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tanp = 1.9, b = 0.1, |^ < 





0.0 

100.0 120.0 140.0 160.0 

M 2 (GeV) 



Figure 2: Mass spectrum of squarks, sleptons, neutralinos, charginos, gluino, 
and Higgs for the case tan/3 = 1.9, fi < 0, and b = 0.1 as a function of the 
wino mass (M2). The line notation is the same as in fig.l. 
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for squarks and sleptons and 

Q x = -2 (18) 

for Higgs doublets. Thus squarks and sleptons obtain a universal positive 
mass squared correction. As a result we can now increase M 2 without the 
light stop mass becoming too small. Then increasing M 2 causes the chargino 
mass to increase. Finally tan/5 can now be increased which raises the 
light Higgs mass above the experimental lower bound. 

In figs. 3(a,b) we plot the sparticle spectrum for tan/5 = 10, M 2 /m z = 2 
as a function of d. We see in fig. 3a that 98 ^ m h < 103 GeV for 0.42 < d ^ 4.0 
(the lightest Higgs mass is given by the solid (dashed) brown line at one loop 
(tree level). Thus the light Higgs problem is solved with the addition of one 
new parameter, d. The Higgs mass is determined by electroweak symmetry 
breaking and is thus only weakly dependent on d. We also see that the mass 
of the other Higgs states, A, H , H + , (evaluated at tree level) decrease as d 
increases, with the pseudo-scalar A remaining the lightest and for d ^ 3.5, we 
have rriA ^ 150 GeV. Squark and slepton masses (fig. 3b), on the other hand, 
increase with d. The light first and second generation squarks and sleptons 
are right handed; the heavy ones are left handed. The third generation 
squarks and slepton eigenstates are mixtures of left and right handed states 
due to A term mixing. The lightest squark is either a stop or sbottom. 
Note, for d < 0.7 the light stop is lighter than the top. This is because the 
right handed stop mass squared is driven negative by renormalization group 
running as a consequence of the large top yukawa coupling. This is the same 
effect which drives the Higgs mass squared negative. Finally chargino and 
neutralino masses implicitly depend on d through their dependence on fi. For 
small d the lightest neutralino and chargino masses are naturally small. The 
lower bound on d is obtained by requiring that the visible width of the Z be 
consistent with the experimental bound (see next section). Recall, the value 
of the gluino mass depends on b. In the figures we have taken b = 0.1; the 
gluino mass scales as b 2 . 

In figs. 4, 5 we show the sparticle spectrum as a function of tan /3 and 
M.2 with the other parameters fixed. In figs. 4(a,b) we take M2/mz = 2, 
d = (0.8,0.5), b = 0.1 and fi > 0. We see (fig. 4a) that for 4 < tan 0< 39 the 
light Higgs mass bound is satisfied. For low and moderate values of tan/? 
the pseudo-scalar (A) and charged Higgs (H + ) are significantly heavier than 
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250.0 



M 2 /m z = 2, tan(3 = 10, b = 0.1, jx > 
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0.0 1.0 2.0 3.0 4.0 

d 

(a) 

Figure 3 (a): Mass spectrum of the light neutral Higgs h (solid (dashed) 
brown) at one loop level (tree level), pseudoscalar boson (A) (red), heavy 
neutral Higgs (H ) (blue), and charged Higgs (H + ) (black) for the case 
M 2 /mz = 2,tan/9 = 10,6 = 0.1 and fi > as a function of d. 
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M 2 /m z = 2, tanp = 10, b = 0.1, |i > 
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500.0 



400.0 



300.0 



200.0 



100.0 




0.0 ' 1 1 1 — 1 1 1 1 — 1 1 1 1 — 1 1 1 1 — 1 

0.0 1.0 2.0 3.0 4.0 

d 

(b) 

Figure 3 (b): Mass spectrum of up squarks (green), down squarks (red), 
sleptons (pink), sneutrinos (orange), neutralinos (blue), charginos (black), 
and gluino (dot-dashed green) for the case M 2 /m z = 2,tan/3 = 10, b = 0.1 
and /i > as a function of d. The first and second (third) generation squark 
and slepton masses are given by solid (dashed) lines. 
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M/n^ = 2, d = 0.8, b = 0.1, n > 
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tanp 
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Figure 4 (a): Mass spectrum of squarks, sleptons, neutralinos, charginos, 
gluino, and Higgses for the case d = 0.8, M 2 /m^ = 2, b = 0.1, and /x > as 
a function of tan/3. The line notation is the same as in fig.l. 
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M 2 /m z = 2, d = 0.5, b = 0.1, \i > 
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M 2 /m z = 2, d = 0.8, b = 0.1, [i > 
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tanp 

(C) 

Figure 4 (c): The mass spectrum of the heavy states (left-handed states) 
of the squarks and slepton (the line notation is same as in fig.3(b)), and 
peudoscalar boson (solid brown line), heavy neutral Higgs (dotted brown), 
and charged Higgs (dashed brown) for the case M 2 /m^ = 2,d = 0.8, b = 0.1, 
and ji > function of tan (5. 
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rriz- Hence the phenomenology of the light Higgs is identical to that of the 
SM. However for large tan/3 ~ 30 the bottom yukawa coupling is relevant. 
In this case, the pseudo-scalar mass m\ = m 2 Hu + m 2 H d becomes small (fig. 
4c). The upper bound on tan f3 is of order 38, since in this case both and 



rriA are greater than ~ 75 GeV, the present experimental bound|T5|. In 
figs. 5(a,b) we take tan/3 = (10, 30), d = 0.8, b = 0.1 and jj, > 0. We see that 
the light Higgs mass increases slowly as M 2 increases with ^ 115 GeV. 
Squark, slepton and gaugino masses also increase as M 2 increases. This is 
because all SUSY breaking masses scale with A. For low values of M 2 of 
order 100 GeV, the lightest neutralino and chargino masses become small; 
the lower bound on M 2 is again determined by the visible width of the Z. 

The gluino mass is arbitrary. Its value is bound from below by about 
1 GeV since, in addition to the GMSB contribution to its mass (eqn. |5|), 
it receives a dynamical QCD contribution to its mass and also radiative 
corrections at the weak scale dominated by the stop; both of order a GeV. 
The gluino mass also increases as M 2 increases and for very large M 2 of order 
a TeV (fig. 6) the gluino mass is above 100 GeV. 

For a large range of the parameter b, either the gluino or the gravitino 
is the LSP.f] In the case that the gravitino is the LSP, then the gluino is 
unstable with the decay rate for a gluino decaying into a goldstino (G) and 
a gluon given by 



5 / Q \ 3 

m- / ., mi 



G + 9) = Te% (19) 



or 



CT; 



g .-lO^^y^cm (20) 



where in eqn. |20] we have neglected the gravitino mass. Hence for all labo- 



ratory experiments the gluino may be considered stable, decaying outside the 

5 Clearly if b is increased, we could imagine a case in which the gluino is heavier than 
the lightest neutralino which is also somewhat light, particularly for small values of M-i 
of order 100 GeV. In this case the neutralino is the LSP (or NLSP). The gluino will then 
decay into the lightest neutralino. This case is highly constrained by standard missing 
energy searches for gluinos. However if the neutralino-gluino mass difference is only a 
few GeV, then the missing energy signal is suppressed and it may not be ruled out. This 
possibility requires further study. If, in addition, the neutralino is also the LSP (lighter 
than the gravitino) it would be stable and would make an excellent dark matter candidate. 
We do not consider this possibility further in this paper. 
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tanp = 10, d = 0.8, b = 0.1, fi > 
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Figure 5 (a): Mass spectrum of squarks, sleptons, neutralinos, charginos, 
gluino, and Higgses for the case tan (3 = 10, d = 0.8, b = 0.1, and /i > as a 
function of M2. The line notation is the same as in fig. 1. 
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tanp = 30, d = 0.8, b = 0.1, n > 
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Figure 6: Same as fig.5 (a 
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detector. Thus the standard missing energy signature for SUSY is gone. The 
gluino will hadronize; forming a hadronic jet. For a heavy gluino, the jet 
may contain only the lightest stable hadron containing the gluino. In any 
event, some fraction of the gluino kinetic energy will be visible in a hadronic 
calorimeter. 

We assume that the LSP is a gluino-gluon bound state, a glueballino 
(-Ro)-Q It is (for all practical purposes) stable because of a conserved R 
parity. Rq will interact in a hadronic calorimeter with a strong interaction 
cross-section. The dominant process at small momentum transfer is governed 
by Regge exchange. In a recent paper by Baer et al. |T7| an estimate for the 
energy loss of Rq in a hadronic calorimeter was obtained using the triple 
pomeron amplitude for the single diffractive process Rq N — > Rq X where 
X denotes the inclusive sum over all final states. It was found that a 25 
GeV gluino with (3 ~ 1 would deposit less than 20% of its kinetic energy 
in the detector assuming 8 hadronic collisions. An Rq, however, could also 
have a significant charge exchange cross-section, given by a triple- Reggeon 
amplitude for the single diffractive process Rq N — > p + X. In this case 
there would be an additional ionization energy loss in the detector, assuming 
p + is sufficiently long lived. Baer et al. have assumed that at each hadronic 
collision the light brown muck surrounding the gluino is stripped off and then 
the gluino re-hadronizes. They parametrize the probability for the resulting 
gluino bound state to be a charged p + by P. We shall discuss their results 
in the next section. It suffices now to say that their results are sensitive to 
this universal parameter P. Clearly a better theoretical estimate of Rq N 
scattering is needed. It is also clear that detailed detector simulations are 
necessary in order to confirm their results. 

We can now test the theory with regards to LEP, CLEO and Tevatron 
data. Note that flavor changing neutral current processes are naturally sup- 
pressed within our framework of combined gauge-mediated and D term SUSY 
breaking. 



6 The other likely candidate LSP is a gluino- (u m — dd)s bound state, the neutral com- 
ponent of an isotriplet p. We shall not consider this possibility further in this paper, see 
for eg. H . 
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3.3 Experimental tests 

First let's consider some tests from LEP. Since the LSP in this theory is a 
gluino, charginos and neutralinos can now decay into qqg. Moreover since 
charginos and neutralinos are relatively light they can be produced at rela- 
tively low center of momentum energies at LEP. We have considered lim- 
its from two processes - - the visible width of the Z (Ar^ (visible)) and 
e + e~ — > hadrons. In figs. 3-6, the lower bound on (d, tan/?, M 2 ), respec- 
tively is determined by the allowed contribution to ATz (visible) < 23.1 x 10 -3 
GeV[|18j. This provides a weak lower bound on neutralino and chargino 



masses. A better bound can be obtained by e + e~ — > hadrons. In figs. 
7(a,b) we give the constraints on the chargino mass coming from the pro- 
cess e + e~ — > hadrons.[] The horizontal line is the OPAL 2 a bound on 
new physics at y/s = 172 GeV[19|]. In fig. 7(a) we show the contribution 



to e + e~ — > hadrons due to chargino and neutralino production times the 
branching fraction for quark decay as a function of the chargino mass (vary- 
ing d) for fixed tan/5 = 10, M-ijmz = 2, fi > 0, b = 0.1. We find a lower 
bound d > 0.44. This translates into a lower bound on the lightest chargino 
mass m~+ > 66 GeV (fig. 7(a)) and neutralino mass m^o > 34 GeV (fig. 3). 
In fig. 7(b) we show the same constraints now varying M 2 instead of d. For 
d = 0.8 and M 2 > 122 GeV, we find m-+ > 71 GeV and > 27 GeV. 

An additional test from LEP is on the number of 4 jet events. The 
dominant decay mode for down squarks is d — > dg. In fig. 4(b) we see that 
the bottom squark can be as light as 75 GeV for large tan/3. Nevertheless 
the calculated cross section for e + e~ — > bb is significantly below the OPAL 
bound for excess four jet events at ^Js = 184 GeV|j20"|. Four jets are also 



possible via the direct process e + e — > qqgg. Baer et al.[fT/f| find strong 



limits using OPAL and L3 data|21] for e + e~ — > Z with Z — > (2 or 4) jets + 
missing energy. This data was used to place limits on the CMSSM process 
e + e~ — > Xi X2 "with the subsequent decay x\ — > q q Xv They argue that a 
gluino mass in the range 5 < m- g < 25 GeV is ruled out. This result seems to 
be relatively insensitive to the probability P. Results from LEP2 at higher 
center of mass energies do not improve upon this limit. 

We have evaluated the rate for b — > 57. The total contribution for fi < 

7 Here we assume that the gluino produces a visible jet. 
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Figure 7 (a): The cross section of e + e~ — *> hadrons at center of mass energy 
of 172 GeV from the chargino and neutralino pair production for the case 
M 2 /mz = 2, tan /? = 10, 6 = 0.1, and \x > as a function of the chargino mass 
(varying d). We also show the OPAL 2 a bound on new physics (dot-dashed 
line) . 
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tanp = 10, d = 0.8, b = 0.1, |x > 
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Figure 7 (b): Same as fig. [7 (a)| for the case tan/5 = 10, d — 0.8, b = 0.1, and 
fi > as a function of the chargino mass (varying M2). 
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M 2 /m z = 2, tan(3 = 10, b = 0.1, \l > 
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Figure 8 (a): The ratio of b — > amplitude of SUSY contribution to Stan- 
dard Model one. The colored lines represent the contribution of charged 
Higgs-top (red), chargino-stop (blue dashed), gluino-sbottom (pink long 
dashed), neutralino-sbottom (green dot-dashed), and total (all SUSY con- 
tribution plus SM contribution) (black solid line) for M 2 /m^ = 2,tan/3 = 
10, b = 0.1, and /i > as a function of d. 
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(b) 

Figure 8 (b): Same as (a) for M 2 /m z = 2,d = 0.8, b = 0.1, and ji > as a 
function of tan f3. 
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Figure 8 (c): Same as (a) for tan/5 = 10, d — 0.8,6 = 0.1, and \i > as a 
function of M 2 . 
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tanp = 30, d = 0.8, b = 0.1, [i > 




(d) 

Figure 8 (d): Same as (a) for tan/5 = 30, d — 0.8,6 = 0.1, and /j > as a 
function of M 2 . 
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is always above the SM result, approaching it asymptotically for large M2.0 
The results for \i > are given in figs. 8(a,b,c,d). It is clear that there 
is a large range of parameters for which the total contribution is below the 
SM value. Our result however only includes one loop analysis and should 
be taken only as an indication that it is possible to obtain results consistent 
with the data fH . 

For d < 0.7 (fig. 3) or tan/3 < 3.8 (fig. 4a) or for any value of tan/? and 
d = 0.5 (fig. 4b), we have rh tl < m t . Then for a sufficiently light gluino, the 
top can decay into a stop + gluino. This range of parameters is testable at 
the Tevatron. In the case that the stop is heavier than the lightest chargino, 
the dominant decay mode of the stop is t\ — > bxt ■ In this case stop decay 
can mimic a top signal. We are now studying this situation. 

It has also been argued by Baer et al.[]i~7| that CDF data for jets + 
missing energy places stringent limits on the gluino mass. These results 
however are sensitive to the parameter P. For P < 1/2 they rule out a 
gluino with mass in the range 20 ^m^^ 140 GeV. However for P = 3/4 
there is an allowed window for the gluino mass 25 < m- g < 35 GeV. Clearly 
these results are more sensitive to the details of how an Rq or p + interacts in 



the detector. For example, as discussed in [I7|], a p + has some probability of 
being mis-identified as a muon if it is seen in the muon tracking chambers. It 
is important that this analysis be redone by both CDF and DO. In addition, 
Baer et al. have explicitly assumed that squarks are very heavy and hence do 
not contribute to the jet + missing energy signal. Inclusion of squarks can, 
via coherent interference, decrease the gluino production cross section and 
hence the rate for jets + missing energy. At the same time, squark production 
and subsequent decay can increase the jets + missing energy rate. Hence the 
Baer et al. analysis must be redone for the model presented in this paper. 



4 Conclusions 

In this paper we have presented a phenomenologically acceptable model with 
a gluino or gravitino LSP. In case the gravitino is the LSP, then the gluino 
is the NLSP. It can decay but only on a cosmological time scale. Thus in 
either case the gluino is stable with respect to any accelerator experiment 
searching for SUSY. We have also discussed some experimental bounds and 

8 This is the reason hgs. 3-6 are given only for /i > 0. 
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possible tests of the theory at LEP, CLEO and the Tevatron. We believe 
that further experimental studies are now justified. 

The theory includes a combination of gauge mediated and D term SUSY 
breaking. The latter contribution was shown to be necessary in order to 
obtain a Higgs with mass consistent with present experimental bounds. 

The gluino is naturally light in this theory because the 10 of Higgs (in 
SO(10)) mixes with the 10 of messengers. It is light because of two effects. 
Firstly, the SU(3) triplet messengers are heavier than the electroweak dou- 
blet messengers; suppressing SUSY breaking mediated by color interactions. 
Secondly, R symmetry (which keeps gauginos massless) is broken at the elec- 
troweak doublet messenger scale; providing an additional suppression. 

There are several other consequences of a theory with a gluino LSP which 
may have cosmological and/or astrophysical significance. It has been shown 
that electroweak baryogenesis is not possible in the SM. On the other hand, 
the necessary conditions for electroweak baryogenesis in the MSSM requires 
irih ~ 100 GeV and < 0[^3|. These conditions are naturally satisfied in 
our model. They gaurantee that there is a first order phase transition which is 
sufficiently strong to shut off baryon violating interactions inside the bubbles 
of broken electroweak phase. If there is also sufficient CP violation; then 
baryogenesis is possible. 

In addition, the Ro provides a natural particle physics candidate for the 
source of ultra high energy cosmic rays, the UHECRON.[] It has recently 
been argued [24| that an Rq with mass in the range from 2-50 GeV could 
reproduce the highest energy cosmic ray shower observed by the Fly's Eye 
Detector . 



Finally, the Ro LSP is not a dark matter candidate since the annihilation 
cross section is too large and hence the relative abundance of Rq to baryons 
is too small.0 There are nevertheless severe constraints on the abundance 
of RqS coming from searches for anomalous heavy isotopes of hydrogen and 
oxygen [2"E|. If R binds to hydrogen, then any R with mass greater than 2 
GeV is ruled out. In addition, if the concentration of R in oxygen is greater 



9 This is even possible if the gluino is the NLSP and decays into a gravitino LSP. For 
example, a gluino with mass ~ 25 GeV and energy characteristic of the highest energy 
cosmic rays ~ 10 8 TeV can travel ~ 10 5 Mpc before it decays. 

10 The uncertainty in the annihilation cross section is quite large ~ 10 . Nevertheless 
the relative abundance of gluinos to baryons is of order (10 -10 — 10~ 7 )(mg/lGeV). For a 
recent re-analysis, sec Bacr ct al.flTfl 
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than one part in 10 16 — 10 19 , it is also ruled out. Note however, if the gravitino 
is the LSP and Rq is the NLSP, these limits are evaded. 
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